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Abstract 

A method is suggested for the computation of the generalized dimensions of fractal 
C^l ■ attractors at the period-doubling transition to chaos. The approach is based on an 

eigenvalue problem formulated in terms of functional equations, with a coefficient 
expressed in terms of Feigenbaum's universal fixed-point function. The accuracy of the 
results is determined only by precision of the representation of the universal function. 



PACS numbers: 05.45.Df, 05.45.-a, 05.10.Cc 

The multifractal or thermodynamic formalism is an important tool for description of 
' strange sets arising in dynamical systems in different contexts. Its basic ideas have been 

0^ ■ clearly formulated e.g. in the paper of Halsey et al. |D. Some of the examples presented 

by these and other authors relate to the fractal attractors that occur at the onset of chaos 
via period doubling and quasiperiodicity The multifractal analysis reveals global 

scaling properties of these attractors, such as the generalized dimensions and the /(a) 
spectra. They are of principal interest because of their universality for systems of different 
nature. Moreover, they allow a measurement in physical experiments J?||. 

One of the well-studied multifractal objects is the Feigenbaum attractor, which oc- 
curs at the period-doubling transition to chaos in unimodal one-dimensional maps with 
quadratic extremum and in a wide class of more general nonlinear dissipative systems 
||, [|. Beside the original procedure of Halsey et al. (namely the construction and analysis 
of the partition functions defined as sums over some natural covering of the attractor), 
several other approaches to the computation of the multifractal characteristics have been 
developed. Bensimon et al. |J used a method based on a break up of a partition sum into 
two components with subsequent use of the scaling property. Kovacs Q suggested a pro- 
cedure of extracting the dimensions from the eigenvalue problem for the Frobenius- Perron 
operator. Christiansen et al. ||] exploited the idea of approximating the strange sets by 
periodic orbits and expressed the desired quantities in terms of cycle expansions. (To 
our knowledge, the calculation of the Hausdorff dimension of the Feigenbaum attractor in 
Ref. || remains the most precise to date.) 

In some sense, the global description of scaling properties in the multifractal formalism 
seems opposite to the local description in terms of the Feigenbaum renormalization group 
approach ||]. The latter is based on the solution of the functional fixed-point equation 
and associated with scaling relations for the evolution operators in a neighborhood of the 
extremum of the map under consideration. 
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In this note we present a novel method for precise computation of the multifractal 
characteristics: the problem may be presented in terms of the Feigenbaum renormaliza- 
tion transformation applied to some auxiliary function. The desired quantities, such as 
the generalized dimensions and the f(a) spectrum, can be extracted from an eigenvalue 
problem, formulated as a functional equation involving Feigenbaum's universal fixed-point 
function. An analogous approach was previously suggested in application to the problem 
of the effect of noise onto the period-doubling transition JO], ^] . That problem appears to 



be linked with one special generalized dimension, as noted e.g. in Refs. J2|, 11 1, and this 
circumstance obviously supports a possibility of the generalization we undertake here. A 
similar idea was used in Ref. 0] for a study of scaling regularities in the Fourier spectrum 
and response function of the quadratic map at the period-doubling transition to chaos. 



Using the representation of the Feigenbaum function from Ref. [13| we obtain the 
generalized dimensions with high precision in excellent agreement with the previously 
known data. Also we present accurate results for generalized dimensions of the multifractal 
attractors at the onset of chaos in the unimodal maps of degrees 4, 6, and 8, for which the 
universal functions g(x) are available in the literature in the form of numerically found 
polynomial expansions fi~4| . 

To estimate the multifractal characteristics for the Feigenbaum attractor by the stan- 

dard approach the generalized partition functions Th(q, r) = Pi/lJ are exploited. Here 

i=l 

q and r are some real parameters, pi = 2 _fc , = \xi — x i+2 k\, and the sequence Xi 
results from iterations of the unimodal map at the limit point of the period-doubling ac- 
cumulation, starting from the extremum point. Obviously, Tk(q, r) = 2~ qk S 2 k(r), where 

S 2 ic(t) = ^ T ■ For each given q an appropriate value of r = r(g) may be found that 

i=l 

ensures an asymptotic equality Tk+i(q, r) = Tk(q, r) as k — > oo. Vice versa, for a given r 
we can select a respective value of q = q(r). This relation of q and r is used then to obtain 
the generalized dimensions and /(a) spectrum. 

For large k the lengths of the intervals li are small, and they may be estimated via the 
derivatives as 

h^ldxi/dxilh. (1) 

In this approximation we can compute them together with the sums S step by step via 
simultaneous iterations of the relations 

Xi+l = f(xi), 

h +1 = If'ixMk, (2) 
S i+ i = S i + l^(x i ). 

We have introduced here an auxiliary function ^(x) which at this moment is supposed to 
be identically equal to 1. 

By twofold application of Eqs. (Q) we obtain 

Xi+2 = f(f(xi)), 

h +2 = \f(f(x t ))f'(x l )\k, (3) 
s t+2 = s l + 1~; 2 [|/'(/(si))r*(si) + *(/(*<))] • 

Now we perform Feigenbaum's scale change x i— > x/a, I i— > l/\a\ (a is the Feigenbaum 
constant) and arrive at the equations, which are of the same form as (|2|), but with new 
functions f ne w(x) = af(f(x/a)), x I / new (x) = Lf )T ^/(x) where Lj jT is the linear operator 

L f>r : ^ \a\ T [\f'(f(x/aW*(x/a) + *(/(x/a))] . (4) 
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This transformation may be repeated again and again. Asymptotically, f(x) con- 
verges to the fixed-point function satisfying the Feigenbaum-Cvitanovic equation g(x) = 
ag(g(x/a)), and the sequence ^fc(x) will follow the recursion ^k+i(x) = Lg^^x). 

Thus, as k — > oo, the function ^ tends to the eigenfunction associated with the largest 
eigenvalue of the linear operator L gT , the eigenproblem being 

u(t)*(x) = \a\ T [\g'(g(x/a))\ T ^(x/a) + V(g(x/a))} . (5) 

(A particular case of this equation for r = 2 appears in the theory of effect of noise onto 
the period-doubling transition. A possibility of computation of the noise scaling constant 



via sums of the derivatives over the Feigenbaum attractor was noted e.g. in Refs. [15, 11].) 

^From the construction, we see that the eigenvalue u(t) indicates a rate of growth or 
decrease of the sums S: 

S 2k (r)^u(r) k . (6) 
To have — ► const as k — » oo, we must have 

v(j) = 2 9 , or q = \og 2 v(r). (7) 

Then, in accordance with the multifractal formalism, we can obtain the generalized di- 
mensions 

and the /(a) spectrum as an implicitly defined relation between the variables 

a = ^— and f = q— r. (9) 

dq dq w 

Although our argumentation starts from the approximate relation ([!]), we believe that 
the final Eq. (||) is exact. The data from the numerical computations presented below 
strongly supports this conjecture: the generalized dimensions are in excellent agreement 
with the best known numerical results, up to all reliable digits. Apparently, in the asymp- 
totics of k — > oo the approximate nature of (|l]) becomes inessential. One might hope that 
a rigorous proof can be found. 

We have performed numerical solutions of the eigenvalue problem (||) for the classic 
Feigenbaum attractor of the quadratic map and for unimodal maps of even integer degrees 
d = 4, 6, and 8. In principle, the achievable precision of the results is determined only by 
accuracy of the approximation of the universal functions. 

With the known polynomial approximation of g(x) and value of the scaling constant a 
we have numerically performed the functional transformation defined by the right-hand 
side of Eq. @. The unknown function ^(x) is represented by a table of its values at 
the nodes of a one-dimension grid on the interval [0, 1] and by an interpolation scheme 
between the nodes. In actual computations it was convenient to use a grid of constant 
step along the axis of variable y = \x\ d and a fourth-order interpolation in terms of y. 
Given the input table for ^(x) the program yields an analogous table as output. 

Suppose we fix r and wish to estimate q. We define an initial condition as ^{x) = 1, 
perform the functional transformation, and normalize the resulting function as ty°(x) = 
^(x)/fy(0). The new function is taken as the initial condition and so on. This operation 
is repeated many times, until the form of the function \&(x) stabilizes. Then, the value of 
^(0) before the normalization is taken to be v{r), and we finally set q{r) = log 2 v(r). 
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To find r for a given q we use the above procedure together with a simple iteration 
scheme for the numerical solution of the algebraic equation q(r) = q. We may then 
calculate D q = r/(q — 1) at q ^ 1. In particular, Dq is the HausdorfF dimension, and D2 
is the correlation dimension. 

To obtain the information dimension D\ it is necessary to determine the limit as q — ► 1, 
that is, as r — > 0. Formally, this follows from L'Hopital's rule: D\ = lim = \^r) = 

g— >1 1 1 \ a 1/q=l 

(<fr) • ^° com P u te this without loss of accuracy we use the following algorithm. For 

t <C 1 let us write ^fc(x) = 2 k \a\ kT [l+Thk(x)] and substitute this expression into Eq. (??). 
To first order we have 

h k +i(x) = \ [h k {x/a) + h k {g{x/a))] + \ In |</( 5 (x/a))| . (10) 

Numerically, representing hk{x) by a table of its values and performing a large number of 
steps of the transformation one can observe that hk+i(x) — hkix) — > = const as k — > 00. 
This implies that * fc oc |a| fcr 2 fc e fceT . On the other hand, ^ fc oc 2 fc "M 2^+^^ = 
2*(i+t/D0. Hence, 

For quadratic maps we have performed the computations based on a polynomial rep- 



resentation of g(x) with coefficients taken from the paper of Lanford [13|. His data are of 
very high precision, but in our calculations the accuracy is limited due to a use of standard 
double-precision arithmetic. As a result, we get not more than 14 true digits in the gener- 
alized dimensions. These data are presented in the first column of Table 1. Note excellent 
agreement of the Hausdorff dimension (up to the last decimal digit!) with the result of 
Christiansen et al. ||. Other dimensions for the Feigenbaum attractor were presented by 
Kovacs pi, and they coincide with our results up to the 10-th digit, the accuracy achieved 
in that work. 

In the remaining three columns of the Table 1 we present results for the generalized 
dimensions of the multifractal attractors at the onset of chaos in unimodal maps of degrees 
4, 6, and 8 obtained using the universal functions given e.g. in Ref . [T^j . 

As an alternative to the traditional definition of the generalized dimensions D q one 
might consider a family of dimensions enumerated by the index r. Let us designate them 
as D T : D T = D q r T \ = r/(g(r) — 1). As mentioned, for r = 2 the equation 
a form studied in the theory of noise effect onto the period-doubling transition 
noise scaling constant is defined as 7 = y/v(2). Hence, the dimension D 2 is linked with the 
effect of noise. The scaling factor 7 is expressed via the dimension D 2 as 7 = 2 1 / D + 1 / 2 . In 
Table 2 we present high-precision data for values of q(2), dimensions D 2 , and noise scaling 
factors obtained from the numerical solution of the eigenproblem (||) for maps of degree 2, 
4, 6, and 8. The data for the factors 7 improve the previously known results [jlq] . (Observe 
that value of q depends on the degree, so one cannot speak on a definite dimension from 
the family D q associated with the noise effect!) 

For the problem of unimodal maps, our method of calculation of the generalized di- 
mensions does not have obvious computational advantages over those of Refs. [| |, |, 
but it does represent the problem in a new light and indicates novel links between global 
and local descriptions of the scaling regularities. An analogous use of renormalization- 
group equations will be feasible in the calculation of multifractal properties in many other 




situations at the onset of chaos, e.g. in bimodal one-dimensional maps [17], asymmetric 
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one-dimensional maps p8| , two-dimensional period-doubling maps JT^], quasiperiodically 



forced maps pQ], and complex analytic maps [21|. Such an approach will be useful es- 
pecially for situations where computations based on the traditional partition-function 
approach are difficult. 
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Table 1: Generalized dimensions for fractal attractors at the onset of chaos in unimodal 
maps of degree d 





d = 2 


d = 4 


d = 6 


d = 8 


D 5 


0.45392270234470 


0.407695571 


0.373232166 


0.351475400 


D 4 


0.46615155691823 


0.426832904 


0.392400635 


0.370142909 


D 3 


0.48077684940009 


0.454569793 


0.421638052 


0.399231039 


D 2 


0.49783645928917 


0.495316676 


0.468035066 


0.447019466 


D x 


0.51709757255124 


0.555181822 


0.544847134 


0.531111008 


Do 


0.53804514358055 


0.642575065 


0.683433256 


0.707102082 


D-i 


0.55991291016494 


0.763919555 


0.946229117 


1.146118382 


D-2 


0.58173600034603 


0.894257449 


1.205507002 


1.510079742 


D-3 


0.60247817187829 


0.992066238 


1.354808070 


1.698747772 


D-4 


0.62126594260209 


1.056616863 


1.445090859 


1.811996998 


D-5 


0.63760518368338 


1.100453275 


1.505301852 


1.887496869 



Table 2: Generalized dimensions D 2 and noise scaling factors for unimodal maps of degree 
d 





d = 2 


d = 4 


d = 6 


d = 8 


9(2) 


5.45324245756108 


6.086657808 


6.654767241 


7.070578662 


D' Z = D q(2) 


0.44911096107158 


0.393185481 


0.353683877 


0.329457884 


7 


6.61903651081803 


8.243910853 


10.037886410 


11.59386214 
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